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Summary: - For the one-dimensional Helmholtz equation we write the corresponding 
time-dependent Helmholtz Hamiltonian in order to study it as an Ermakov problem and 
derive geometrical angles and phases in this context. 

Our starting point is the Helmholtz equation in the form written by Goyal et al Q 
and by Delgado et al 

+ X^{x)i,{x) =0, (1) 

i.e., as a Sturm-Liouville problem for the set of eigenvalues A G i? defining the 
Helmholtz spectrum on a given, closed interval [a,b] of the real axis on which a non 
trivial V' is vanishing at both ends (Dirichlet boundary conditions). Eq. (1) occurs 
for example in the case of TE (transverse electric) mode propagation in planar 
waveguides with a refractive index which varies continuously in the x direction but 
is independent of y and z and in equivalent problems in acoustics. We perform 
the mapping of Eq. (1) to the canonical equations for a classical point particle in 
the known way (for recent discussions see [|[ ||). Let ip{x) be any real solution of 
Eq. (1). Defining x=t, — q, and — p, Eq. (1) becomes 

S ^ ^ (^) 

f = (3) 

with the spectral condition q{a) — q(b) — 0. Eqs. (2,3) are the canonical equations of 
motion for a classical point particle as derived from the time-dependent Hamiltonian 

H{t) = ^ + Xq^itf^ . (4) 

2 

For this Hamiltonian the triplet of phase-space functions Ti — T2 — pq, and 

2 

T3 = ^ forms a dynamical Lie algebra (i.e., H = J2n hn{t)Tn{p, q)) which is closed 
with respect to the Poisson bracket, or more exactly {Ti,T2} = —2Ti, {T2,T3} ~ 
— 2T3, {TijT^} = —T2. The Helmholtz Hamiltonian can be written down a.s H = 
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Ti + X(l){t)T^. The so-called Ermakov invariant / Q belongs to the dynamical 
algebra 

I = Y,^ir{t)Tr , (5) 



and by means of 



one is led to the following equations for the unknown ^r{t) 



fir 



E 



(6) 



(7) 



where C^^n the structure constants of the Lie algebra that have been already 
given above. Thus, we get 



Ai = -2^2 

/is = 2X(l){t)n2 ■ 



(8) 



The solution of this system can be readily obtained by setting /ii — givirLg 
IJ.2 — —PP and /i3 = + where p is the solution of the Milne-Pinney equation 13 



p + \(t){t)p = -3 , 

In terms of the function p[t) the Ermakov invariant can be written as follows 



(9) 



{pp - pqf 



2p2 



(10) 



Next, we calculate the time-dependent generating function allowing one to pass to 
new canonical variables for which / is chosen to be the new "momentum" 



S{q,I,p{t))= dq p{q ,I,p{t)) 



leading to 



9 . 

S{q, I, p{t)) — h /arctan 

2 p 



2p2 



where we have put to zero the constant of integration. Thus 



dS 

a7 



arctan ( 



Moreover, the canonical variables are now 

q = pV2Ismd , 



(11) 



(12) 



(13) 



(14) 
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and 



p — — — ^ cos + pp sin 6 



The dynamical angle will be 

whereas the geometrical angle reads 

= - 
2 



p2 d (p 



2 dt' \p 



dt 



(PP) 



dt 



(15) 



(16) 



(17) 



For periodic parameters il{t), with all the components of the same period T, the 
geometric angle is known as the (nonadiabatic) Hannay angle Q and in terms of p 
is given by 

Aefj ^ ^ (j) pdp . (18) 

Passing now to the quantum Ermakov problem we turn q and p into quantum- 
mechanical operators, q and p — —ih-^, but keeping the auxiliary function p as a 
c number. The Ermakov invariant is a constant Hcrmitian operator 



(19) 



in^iv(g,i)> = ^if + m)f)\mt)) 



and the time-dependent Schrodinger equation for the Helmholtz Hamiltonian reads 

(20) 

The goal now is to find the constant eigenvalues of / 

/|V'n(<?,i) =K„|V„((7,i)) (21) 

and to write the explicit superposition form of the general solution of Eq. (20) |^ 

V'(g,i) = E^»^""^*^^"('?'*) (22) 

n 

where C„ are superposition constants, ipn are the (orthonormal) eigenfunctions of 
/, and the phases a„(t), known as the Lewis phases in the literature j|, 10 , are to 
be found from the equation 



(23) 



The key point for the quantum Ermakov problem is to perform a unitary transfor- 
mation in order to obtain a transformed eigenvalue problem for the new Ermakov 
invariant / = UIU"^ possessing time-independent eigenvalues. It is easy to get the 
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required unitary transformation as U — exp[— |-^^] and the new Ermakov invari- 

/ 2-2 "2 t)2 — 

ant will be / = ^ — h Therefore, its eigenfunctions are cx Hn{9 / yfi), 
where Hn are Hermite polynomials, = and the eigenvalues n{n + \). 

Thus, one can write the eigenfunctions V'n as follows 



V.„oc4exp(ijV)exp( 



(24) 



The factor l/p^^^ has been introduced for normalization reasons. Using these func- 
tions and performing simple calculations one is lead to the geometrical phase 





{pp) - (? 


'to 





dt 



The cyclic (nonadiabatic) Berry phase reads 



n = in+l)<f Pdp 

J c 



(25) 



(26) 



Our last issue is the adiabatic approximation. We have found that a good adiabatic- 
ity parameter is the inverse square root of the Helmholtz eigenvalue, leading to 
a "slow time" variable r — -^t- The adiabatic approximation has been thoroughly 
studied by Lewis . If the Helmholtz Hamiltonian is written as 



(27) 



then the corresponding Milne-Pinney equation reads 

1- ,/N 1 
-p + 4>{t)p= -3 , 

A p-^ 

while the Ermakov invariant will be a 1/VA-dependent function 

{pp - n/^f 



/(1/VA) = 



r 

2p2 



(28) 



(29) 



In the adiabatic approximation, Lewis |^ obtained the general solution of the Milne- 
Pinney equation in terms of the linearly independent solutions / and g of the equa- 
tion of motion \q + Vi^{t)q for the classical oscillator (see Eq. (45) in ^). Among 
the explicit examples given by Lewis, it is Vl{t) — bt"^/^, m 7^ —2, h = const, which 
is directly related to the dielectric planar waveguides since it corresponds to power- 
law index profiles {n{x) cx x"^/"^). For this case, Lewis obtained a simple formula 
for the p function when required to be 0(1) in l/VX 



7i 



(m-K2) 



{V)iif 



(30) 
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where H^p^ and H^p^ are Hankel functions of order /3 = l/(m + 2), y = 
and 71 = ±1, 72 = ±1. Even more useful for the technological application may be 
Lewis's specialization to m = —2^^, n = ±1,±2, ... leading to 



where 



G(i,l/\/A) 



k (n 



ky. 



i/Vx 



k=0 



fc!(n - fc)! V2i6(2n+ 1) 



t (2"+l) 



(31) 



(32) 



Thus, one gets p as a polynomial in the square of the adiabaticity parameter, i.e. 
A~^, which has an infinite radius of convergence. One can calculate the topological 
quantities by plugging the above explicit Milne-Pinney functions into the corre- 
sponding formulas. However, Lewis carried out and checked a recursive formula 
in 1/A to order l/A'^, which can be used for any index profile. It is 



P = Po + Pi/A + P2/A^ + ps/A^ + 



(33) 



with po = rj~^/^ = 0^^/^(a;); for the other coefficients pi see the appendix in ||^. 
The main contribution to the topological quantities comes from po- In the case of 
power-law index profiles one obtains the geometrical angle 



A03 = 



46(m -I- 2) 



t-(f+i) 



-(f +1) 



(34) 



and a similar formula for the geometrical quantum phase. For periodic index pro- 
files, one can work out Hannay angle and Berry's phase according to their cycle 
integrals. Finally, we notice that choosing 4>{x) = <I>(x) -I- ^""^^ , which means non- 
linear waveguides, leads to a somewhat more general time-dependent Hamiltonian 
that has been discussed in the Ermakov perspective by Maamache [|o| . 
In summary, we presented an application of the Ermakov procedure to one-dimensional 
Helmholtz problems in a way which may be of direct experimental relevance. We 
notice that very recently Petrov ||ll| discussed the evolution of Berry's phase in a 
graded-index medium, whereas several years ago, Goncharenko et al studied 
the Ermakov approach in nonlinear optics of elliptic Gaussian beams. 
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